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ABSTRACT 

We performed single point [C I] 3 Pi- 3 Po and CO J=4-3 observations to¬ 
ward three T Tauri stars, DM Tau, LkCa 15, and TW Hya, using the Atacama 
Large Millimeter/submillimeter Array (ALMA) Band 8 qualification model re¬ 
ceiver installed on the Atacama Submillimeter Telescope Experiment (ASTE). 

Two protostars in the Taurus L1551 region, L1551 IRS 5 and HL Tau, were also 
observed. We successfully detected [C I] emission from the protoplanetary disk 
around DM Tau as well as the protostcllar targets. The spectral profile of the 
[C I] emission from the protoplanetary disk is marginally single-peaked, suggest¬ 
ing that atomic carbon (C) extends toward the outermost disk. The detected 
[C I] emission is optically thin and the column densities of C are estimated to 
be < 10 16 cm -2 and ~ 10 17 cm -2 for the T Tauri star targets and the proto¬ 
stars, respectively. We found a clear difference in the total mass ratio of C to 
dust, M(C)/M(dust), between the T Tauri stars and protostcllar targets; the 
M(C)/M(dust) ratio of the T Tauri stars is one order of magnitude smaller than 
that of the protostars. The decrease of the estimated M(C)/M(dust) ratios for 
the disk sources is consistent with a theoretical prediction that the atomic C can 
survive only in the near surface layer of the disk and C + /C/CO transition occurs 
deeper into the disk midplane. 
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Subject headings: circumstellar matter — protoplanetary disks — stars: pre-main 
sequence — submillimeter: stars 


Introduction 


The gaseous component of protoplanetary disks is crucial for understanding the planet 
formation process because it affects the structures and evolution of the disk via chemical 
reactions. Theoretical studies have predicted that protoplanetary disks have multiple layers 
in their vertical direction due to stellar radiation. These include a relatively cool midplane, 
in which most of the molecule s are de picted onto grains, a molecular-rich intermediate layer, 
and a hot surface layer (e.g., iBergin et al. 2007]). Because the physical environment differs 
among layers, gas phase abundances also vary significantly. 

Atomic carbon (C) is considered to be abundant in the disk surface where far-ultraviolet 
(FUV) photons drive the energy balance and gas chemistry, i.e., the photon dominated 
region (PDR). Theoretical studies predict that in the PDR, FUV photons with energies 
less than the ionization energy of hydrogen (13.6 eV) dissociate molecular h ydroge n and 
carbon monoxide (CO) and ionize C to yield ionized carbon (C + ) (Tielens & Holfenbach 
19851 ). When significant attenuation is achieved, a thin H/fR transition layer should appear, 
beyond which hydrogen molecules dominate. Because the ionization energy of C is 11.6 eV, 
C+ is dominant at the lower density regions of the PDR. Deeper to the PDR, the carbon- 
ionizing radiation is attenuated, and CO forms via various chemical reactions. As a result, 
for the uniform PDR, C exists in a thin layer sandwiched between the C + and CO layers. 

B ecause phot oevap oration at the disk surface is believed to govern the gas dispersal 
(Hollenbach et al. 1994, 20001), observations of C are crucial for understanding disk dissi¬ 
pation. Moreover, the observations of C are also important in terms of dust evolution in 
the disk because it affects the chemistry as a result_of changes in the UV radiation field 
propagation and thus the intensity of the C line (IJonkheid et ah 2004). 


However, detection of the submillimeter fine structure C lines has not been reported for 
protoplanetary disks thus far. Although a few studies have attempted to detect the emission 
lines of C from protoplanetary disks around intermedia t e-mass stars, only up per limits have 
been obtained (jChapillon et al.lDOlOt IPanic et ahll2010c ICasassus et al.ll2013lh 


In this paper, we present the results of [C I] 3 Pi- 3 Po and 12 CO J— 4-3 observations with 
the Band 8 qualification model receiver mounted on the Atacama Submillimeter Telescope 
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Experiment (ASTEj^J, and report the first detection of [C I] emission from a protoplanetary 
disk. We estimate the optical depth of the emission to derive the column density of C, and 
discuss possible evolutionary trends in the total mass of atomic C. 


2. Observations 


We performed [C I] 3 Pi- 3 Po (492.161 GHz) and CO J =4-3 (461.041 GHz) observations 
toward three classical T Tauri stars (TTSs), DM Tau, LkCa 15, and TW Hya, in November 
2010 using the ASTE. Target information is listed in Table Q3 In addition, we made [C I] 
observations toward two protostars (PSs) in Taurus, L1551 IRS 5 and HL Tau, to check the 
consistency of the observational data. We also used these data as a reference of the [C I] 
intensity during the PS phase of star formation. 


The observations were taken w ith an ALMA Band 8 (400-500 GHz) qualification model 
heterodyne receiver (Band 8 QM; ISatou et ah 20081 ) mounted on the ASTE. The Band 8 
QM uses sideband separating receivers that detect two orthogonal polarizations and down- 
convert the sideband-separated intermediate frequency signals to 4-8 GHz. The half power 
beam width (HPBW) was 17" and the main beam efficiency //mb was estimated to be 45% 
from observations of Jupiter. For the backend, we used MAC, a 1024 channel digital auto¬ 
correlator, which has a band width of 128 MHz and a resolution of 125 kHz, corresponding 
to 78 and 0.076 km s' 1 * , respectively. 


The position switch method was employed and a single spectrum toward the source 
position was obtained for each target. Because our TTS targets are well known to be less 
affected by surrounding clouds, we selected OFF positions near the targets to effectively 
remove atmospheric fluctuations; 28, 15, and 1 arcmin from the stellar positions of DM Tau, 
LkCa 15, and TW Hya, respectively. For the PSs, we employed an OFF position of («j 2 ooo, 
<5j2ooo) = (4 h 32 m 6!0, 17° 4 8' 51"0), which was determined to avoid the CO J =1-0 emission 


of the L1551 cloud (lYoshida et ah| 2010). 


The telescope pointing calibration was performed every 1.5-2 hours by observing O-Cet 
and IRC+10216 in the CO J =4-3 emission line, and the resulting pointing accuracy was 
<3". Any time variation in intensi ty sc ale was checked by observing part of the Orion Horse- 
head nebula (position A; Philipp et ah 2006 ) two or three times daily and was found to be 


1 The 10 m submillimeter telescope operated by Nobeyama Radio Observatory (NRO), a branch of Na¬ 

tional Astronomical Observatory of Japan, in collaboration with the University of Chile, and Japanese 
institutes including University of Tokyo, Nagoya University, Osaka Prefecture University, Ibaraki University, 

and Hokkaido University. 
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less than ±10%. The [C I] data were obtained under good sky conditions ( 7220 GHz 0.05), 

whereas t 2 2 ogh z =0.06-0.10 for the CO observations. The single sideband system noise 
temperatures were typically 2100 K and 3100 K for [C I] and CO, respectively. 

Data reduction and analysis were performed using the Common Astronomy Software 
Applications (CASA) version 3.3.0 software suite, in addition to its ASAF0 modules. It 
should be noted that we reduced only one polarization dataset because of a beam alignment 
problem during the observations. Bad data were removed by the sdflag task, baseline fitting 
and subtraction of the baseline from the spectra were performed with the sdbaseline task, 
and the data sets were combined in each source with the sdaverage task. The final spectra 
were obtained after smoothing with a five pixel boxcar function along the channel axis by 
the sdsmooth task, corresponding to a velocity resolution of ~ 0.4 km s -1 . 


Results 


3.1. [C I] 3 Pi— 3 P 0 and CO J =4—3 spectra 


Figure [T] shows the [C I] 3 Pi- 3 P 0 and CO J =4-3 spectra toward the three T Tauri 
star targets, and the derived spectral line parameters are listed in Table [2l We successfully 
detected [C I] emission with a peak intensity of 2% = 0.12 K toward DM Tau. The line 
profile is single peaked with a central velocity of Vlsr =6.1 km s _1 , which is consistent with 
the systemic velocity of the circum stellar disk as estimated by molecular line observations 


(5.9 km s , e.g., Handa et al. 19951k The emission appeared to be slightly asymmetric with 


2 ATNF Spectral Analysis Package 


Table 1. Target information 


Source 

R.A. 

(J2000) 

Decl. 

(J2000) 

d 

(pc) 

SR 

type 

M* 

(M 0 ) 

Age 

(Myr) 

M(dust) 
(10- 4 Mq) 

Rdisk(CO) 

(AU) 

incl. 

(°) 

./"disk 

refs. 

DM Tau 

04:33:48.72 

+18:10:10.0 

140 

Ml 

0.5 

4.3 

2.4 

890 

34 

0.32 

1,4,5,6 

LkCa 15 

04:39:17.80 

+22:21:03.5 

140 

K5 

1.0 

3-12 

4.8 

905 

52 

0.25 

1,6,7 

TW Hya 

11:01:51.91 

-34:42:17.0 

54 

K7 

0.8 

10 

3.4 

215 

7 

0.15 

2,8 

L1551 IRS 5 

04:31:34.07 

+18:08:04.9 

140 
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1 

1 

HL Tau 

04:31:38.44 

+18:13:57.7 

140 
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a somewhat larger extent to redshifted v elocities. This is similar to the spectra l shap e of the 
CO J =4-3 and other CO transitions, f Guilloteau fc Dutrev 1994 : Thi et al. 2001). Here, 
the emission component visible ne ar ~ 9 km s -1 was not analyzed because it most likely 
originates from the ambient cloud (Handa et ah 19951). For TW Hya and LkCa 15, we did 
not detect any [C I] emission, and obtained only the 3a upper limits shown in Table [2j To 
estimate upper limits to the total integrated intensity of the [C I] emission lines, we assumed 
velocity widths of 3 and 1 km s^ 1 for Lk Ca 15 a nd TW Hva, respectively, the values obtained 
by single dish CO observations (Ivan Zadelhoff et al] 200 li). 


Significant [C I] emission toward the PSs L1551 IRS 5 and HL Tau was detected, as 
shown in Figured The spectral profiles are single-peaked at values close to their systemic 
velocities. The emission intensities are markedly higher than that of DM Tau. 


Emission in the CO J=4-3 line was detected toward DM Tau and TW Hya. The 
peak velocities and line profiles are in good agreement with those o f high-J CO spectra 
( Guilloteau fe Dutrev 1994 ; Thi et al. 2001 : van Zadelhoff et al. 2001 b The double-peaked 
profile of DM Tau lies close to the systemic velocity. After correcting for the beam size 
and efficiency, the integrated intensity of CO emission from TW Hya was comparable to the 
marginal value reported by iKastner et ah (1997) within the uncertainties. 


3.2. Optical depth of [C I] emission and C column density 


We derived the optical depth of the [C I] emission and thus the C column density using 
the equations reported by Oka et al. ([20011 ). The optical depth at the peak velocity, T[c i], 
is given by 

i f i TI/t)mb \ m 

T [C I] = - 111 1 - -- f77 ^-), (1) 


1 /disk[J(T ex ) - J(T bg )] 

where /disk is the beam filling factor of the source, T ex is the excitation temperature, T bg 
is the temperature of the cosmic background radiation, and J is the radiation temperature 
defined by 

hv/k-Q 


J(T) = 


( 2 ) 


exp(/ri///cBT) — 1 

Here h is the Planck constant, v is the observed frequency, and is the Boltzmann constant. 
The factor /disk was derived from the ratio of the solid angle of the disk to that of the 
telescope’s beam, /disk = Ddisk/^A- The solid angle Ddisk was estimated b y ass uming that 
the extent of the [C 11 emission is identical to that of the resolved CO disk (IPietu et al.ll2007l: 
Andrews et al. 2012). The outer radius of the CO disk and its inclination angle, in addition 
to the estimated /disk, are listed in Table [U We employed / disk = 1 for the PSs because 
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they were expected to be embedded in extended envelopes that extend over the beam area. 
For the TTSs, we assumed a constant T ex of 50 K, as suggested by theoretical studies of 
protoplanetary disks around typical T Tauri st ars that predict the C + /C/CO transition 
layer to appear at regions with T gas ~ 50 K (IJonkheid et al.l 120041: iKamp fe Dullemond 


2004 ). Conversely, we assumed T ex = 20 K for the PSs, which is the accep ted value for [C I] 


observations of a PS surrounded by an infalling envelope (iCeccarelli et al.l 1 9981 ). 


Assuming that local thermodynamic equilibrium (LTE) applies, we derived the averaged 
C column density over the source via 


N(C) = 1.98 x 10 

J(n 8 ) 

J(T, X ) 


15 / ^ A di> 


X 


1 - 


?1mb 

-i 


Q(T ex ) exp 


\kT e J 


T IC I] 

1 - exp(—T[ C I]) 


/disk 


(3) 


Here, Q(T ex ) is the ground-state partition function for neutral atomic carbon, 


Q(T ex ) = l + 3exp(- -J^)+5exp(- , (4) 

where E x and E 2 are the energies of the J =1 and 2 levels, respectively, and T[c ij is the 
velocity-averaged optical depth of the [C 1] emission. In this study, however, we used the 
optical depth at the peak velocity instead of the velocity-averaged optical depth, which gives 
an upper limit to the column density. It should be noted that we assumed the optically thin 
condition (tjc i] Cl) for LkCa 15 and TW Hya because no [C I] emission was detected. The 
C column density so estimated are listed in Table [3] The total mass of atomic carbon, M(C), 
was derived by integrating the column density over the solid angle of the source. Therefore, 
M(C) in the telescope beam can be uniquely determined regardless of source distribution. 


Table [3] clearly shows that all the detected [C I] emission lines are optically thin. The 
estimated C column density of the TTSs is typically ~ 10 16 cm -2 or less, whereas that of 
the PS is one order of magnitude higher (~ 10 17 cm -2 ). Accordingly, the M(C) values of 
the T Tauri star targets are one order of magnitude smaller than those of the PSs. These 
results can be attributed to the lack of dense envelopes around the TTSs. 


4. Discussion 

4.1. Origin of the [C I] emission 

The [C I] emission detected from DM Tau is most likely associated with the rotating gas 
disk around the star. The peak velocity and full width half maximum are consistent with 
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those of the disk in molecular lines (e.g., Saito et al.l 119951 ). Because of the modest upper 
energy of the ground state atomic C fine structure line (~ 24 K), the [C I] emission is likely to 
be weighted toward the emission from the outer most part of the disk due to beam dilution. 
In fact, the detected [C I] emission shows a single-peaked profile near the systemic velocity, 
suggesting that the emission may dominate at the region where the rotation velocity is small, 
i.e., the outer part of the disk if Keplerian rotation applies. In addition, the fact that the 
total mass of C is one order of magnitude smaller than those of the PSs (Table |3|) supports 
the disk origin of the [C I] emission. However, the actual nature of the [C I] emission of DM 
Tau remains unclear due to the large HPBW of the telescope. High-resolution observations 
are required to reveal the detailed spatial distribution of the [C I] emission over the disk. 

The estimated [C I] inte nsity of DM Tau is consistent with theoretical models of pro¬ 


toplanetary disk chemistry. Jonkheid et ah ([2004), for example, calculated the [C I] line 


intensity of a T Tauri star with a massive disk of 0.07 M 0 through a self-consistent treat¬ 
ment of the PDR near the disk surface. The resultant peak intensity is expected to be ~0.1 
K in after correcting for the ASTE beam size and the distance to the source, which is 
comparable to our observed intensity. 

The higher intensities of the PSs than those of the TTSs suggest that the atomic carbon 
line originates mainly from the envelopes. Because the velocity width of the [C I] emission 
is much narrower than that of the wing emission of a molecular outfl ow 10 km s^ 1 ), the 


high velocity outflow is unlikely to affect the [C I] emission flMoriartv-Schieven et al.l 12006 


Stoiimirovic et al. 120061) . There is a possibility of a low velocity cavity wall swept up by the 
outflow. High-resolution [C I] observations are required to separate the envelope from such 
a low velocity entrainment component. 


4.2. Difference of M(C)/M(dust) between T Tauri stars and protostars 

We introduce the total mass ratio of atomic C to dust grains, M(C)/M(dust)(= R ), to 
investigate the evolutional variation of M(C). The dust mass, which is generally estimated 
by optically thin (sub-)millimeter continuum observations, reflects the total gas mass of the 
system if the gas-to-dust ratio is constant. Therefore, R is an indicator of a fraction of C 
with respect to the total amount of gas, which is described as R oc g/d x X(C), where g/d 
is the gas-to-dust ratio, and AT(C) is the fractional abundance of C with respect to H 2 . It 
is important to state that, especially for a protoplanetary disk, the atomic C and the dust 
grains trace distinct vertical layers of the disk; the atomic C presumably exists in the disk 
near-surface layers while the dust grains are mainly at the disk midplane. Therefore, R 
practically indicates the mass of C near the PDR with respect to the total mass. 














We compiled M (dust) from previous studies, as listed in Table [TJ All M(dust) data 


except that for TW Hya, were calculated from the total masses in Andrews & Williams 
(120051) . who measured (sub-)millimeter flux densities with single dish telescopes. Their beam 
sizes were typically ~10-15", which is comparable to that of our [C I] observations. These 
authors estimated the total (gas+dust) disk mass by a simple disk model with power-law 
density and temperature profiles under an assumed gas-to-dust ratio of 100. In Table [U we 
show the dust mass after correcting for the assumed gas-to-dust ratio. For TW Hya, we took 
the mass estimate from IWilner et ah (12000!) . who also used a similar simple disk model to 
estimate the disk mass. Although their observations were made with an interferometer, the 
flux density of the entire disk was correctly recovered because the disk size is comparable to 
the detectable scale mentioned in Wilner et ah (120001 ). 


The derived R is listed in Table [31 A clear difference in R is apparent between the 
TTSs and PSs; R is typically an order of ICC 4 or less for the TTSs, whereas the PSs show 
R = 3.4 x 10~ 3 , on average. It should be noted that these results are not significantly affected 
by the difference in the assumed value of T ex . If we use T ex =50 K for the PSs, which is the 
same value as in the TTSs, we obtain R = 2.7 x ICC 3 , on average, a value still significantly 
higher than that for the TTS values. When T ex =10 K, a typical value in molecular cloud 
cores, the difference was more noticeable. 


As defined above, the decrement of R can be interpreted by either a reduced A"(C) or 
g/d. The former case means that 27(C) decreases as star formation progresses. The chemical 
lifetime of C h as be en theor etically predicted to drop substantially with increments in the gas 
density (iLeung et ah 1984 ). Because the gas density significantly increases from a tenuous 
protostcllar envelope to a dense disk near a TTS, a reduction in 27(C) is expected at the 


TTS stage. In fact, the typical density at the disk midp 
at 100 AU if we assume the mini mum mass solar nebu la i 
density of the PSs is 10 5 cm 3 (IMomose et al.lll998l : ISaito et al.1 


ane is estimated to be 


10 9 cm 3 


Havashi 


9811) . whereas the typical 


2001). If the reduced g/d 


is the case, the gas in the disk must be depleted to a low level of g/d < 10 to explain the R 
decrement of at least one order of magnitude. However, such a high depletion of overall gas 
is not expected to occur in the disk, as discussed in ( 14.31 


4.3. Non detection of [C I] emission toward LkCa 15 and TW Hya 

Only upper limits were obtained for the [C I] emission from LkCa 15 and TW Hya, 
despite the fact that they also have disks as massive as that of DM Tau. One possible 
explanation is that overall gas depletion with disk evolution leads to a depletion of C. Since 
all the targets are in the evolutional stage of a so-called transitional disk, their disks are 
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expected to be evolved with respect to a filled disk (Andrews et ah 201 ll 2012). Moreover, 
their advanced stellar ages imply that the disks of LkCa 15 and TW Hya are older than the 
DM Tau disk. The actual gas mass is difficult to directly estimate because of the large optical 
depths of molecular lines. Some modeling efforts by using CO data indicate a lower gas-to- 


dust ratio with respect to the standard interstellar value (IWilliams fe Bestl 12014]; iThi et al. 


2010 ). However, as mentioned in 1 14.21 the gas mass estimation for TW Hya by an HD line 


emission, which is believed to be a good tracer of gas mass, indicates that the disk is massive 


so that the gas-to-dust ratio is comparable to the standard interstellar value (jBer gin et ah 


2013). Taking into account the fact that the disk_ model with a, nominal gas-to-dust ratio of 


100 reproduces well the observed intensity (Jonkheid et ah 2004), an overall gas depletion 
is unlikely the case. An alternative is the depletion of volatiles such as C-containing species 
at the surface with disk evolution. Such a time-dependent sink mechanism of C is proposed 
based on the observational data of HD and CO isotopologue lines (Favre et al. 120131 ). 


We acknowledge the ASTE staff for the operation and maintenance of the observational 
instruments. We would like to thank H. Nomura and D. Ishimoto for fruitful discussions. 
T.T. and M.M. were supported by JSPS KAKENH1 grant No. 23103004. Y.S. was supported 
by the French National Research Agency (grant No. ANR11BS560010 STARFICH). 
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Fig. 1.— [C I] (left) and CO (right) spectra toward the three TTSs, DM Tau (top), LkCa 
15 (middle), and TW Hya (bottom). The vertical red line indi cates the systemic velocity of 
the circumstellar disk I Handa et al. 1995 : Qi et al. 2003 . 2004 1. 
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Fig. 2.— [C I] spectra of the PSs, L 1551 IRS5 (left) and HL Tau (right). The vertical red 
lines indicate the systemic velocity (jSaito et a h 20011). 


Table 2. Parameters of [C I] 3 Pi- 3 P 0 and CO J =4-3 spectra 


Source 


[C I] 

3 Pi- 3 Po 



CO J =4-3 


rp* i 

1 A 

(K) 

Vlsr 

(km s _1 ) 

AV 

(km s -1 ) 

J' T X dv 

(K km s _1 ) 

rp* 1 

1 A 

(K) 

Vlsr 
(km s -1 ) 

At/ 

(km s _1 ) 

J' T X dv 

(K km s _1 ) 

DM Tau 

0.12±0.03 

6.1 

1.4 

0.17±0.05 

0.18±0.03 

6.2 

i.i 

0.20±0.05 

LkCa 15 

<0.15 




<0.18 




TW Hya 

<0.24 




0.62±0.08 

2.8 

0.7 

0-43it0.il 

L1551 IRS 5 

1.77±0.10 

6.3 

1.6 

2.77±0.21 





HL Tau 

2.26±0.08 

6.4 

2.0 

4.55±0.18 






1 3a upper limit for non-detection. 


Note. — Column (1) gives the source name. Columns (2) and (6) give the peak intensity and root-mean-square noise level of the line 
in TJ. Columns (3) and (7) give the local standard-of-rest (LSR) velocity at the emission peak. Columns (4) and (8) give the velocity 
width in full width at half maximum (FWHM) of the line. Columns (5) and (9) give the total integrated intensity of the line and its 
uncertainty. 
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Table 3. Physical parameters derived from [C I] line. 


Source 

yc i] 

N( C) 

(10 16 cm -2 ) 

M(C) 
(10- 7 M 0 ) 

M(C) 

M (dust) 

(io- 4 ) 

DM Tau 

0.02 

1.6±0.5 

0.75=t0.22 

3.1 

LkCa 15 

<0.03 

<4.1 

<1.48 

<3.1 

TW Hya 

<0.09 

<3.6 

<0.12 

<0.4 

L1551 IRS 5 

0.47 

10.7±0.8 

15.3±1.2 

31 

HL Tau 

0.65 

19.0±0.8 

27.4±1.1 

46 


Note. — Column (1) gives the source name. Column (2) gives the 
optical depth of the [C I] line. Column (3) gives the column density 
of C after correction for the beam filling factor of the disk. Column 
(4) gives the total mass of C. Column (5) gives the ratio of the total 
mass of C to the dust mass listed in Table [T] 








